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Abstract 

A prominent example of a topologically ordered system is Kitaev’s 
quantum double model 2?(G) for finite groups G (which in particular 
includes G = Z 2 , the toric code). We will look at these models from the 
point of view of local quantum physics. In particular, we will review how 
in the abelian case, one can do a Doplicher-Haag-Roberts analysis to study 
the different superselection sectors of the model. In this way one finds that 
the charges are in one-to-one correspondence with the representations 
of T>(G), and that they are in fact anyons. Interchanging two of such 
anyons gives a non-trivial phase, not just a possible sign change. The 
case of non-abelian groups G is more complicated. We outline how one 
could use amplimorphisms, that is, morphisms 2t ^ M„(2t) to study the 
superselection structure in that case. Finally, we give a brief overview 
of applications of topologically ordered systems to the held of quantum 
computation. 


1 Introduction 

A fundamental result in local quantum physics is that in high enough dimen¬ 
sional space-times (or, for well enough localised particles), charged particles are 
either (para-)bosons or (para-)fermions |19] . This is no longer true in lower di¬ 
mensional space-times |25] . Instead of a representation of the symmetric group, 
representations of the braid group may be obtained from interchanging identi¬ 
cal particles. Even though our world appears to have three spatial dimensions, 
many systems effectively behave like two or one dimensional systems, opening 
the possibility that they may be described by a low dimensional effective theory, 
with excitations with braided statistics. 

Interest in such systems has sparked in recent years, getting attention from 
the theoretical and mathematical physics communities, condensed matter physi¬ 
cists, quantum information theorists, and mathematicians. Part of the reason 
for this is the discovery of topological quantum computation (first proposed in¬ 
dependently by Kitaev m and Freedman [27]), a field at the intersection of 
quantum theory, quantum computation and quantum topology (61] . One of the 
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attracting features is that the use of topological properties of the system can 
lead to much better fault-tolerance with respect to (local) perturbations. An 
overview and candidates for systems that can be used for topological quantum 
computation can be found in the review m- 

On the mathematical level the behaviour of the anyons, the quasi-particle 
excitations with braided statistics, can be captured by the concept of a braided 
fusion category. This category encodes all information on how two anyons can 
combine (“fusion”), what happens if we interchange them (“braiding”), which 
types can exist in the system, and so on. This structure is in fact very well 
known in local quantum physics. One of the highlights of algebraic quantum 
field theory is the study of superselection sectors initiated by Doplicher, Haag 
and Roberts [IIIID], see also this volume [Ml [55] . This leads in a natural way 
to a fusion category as above. Also in the local quantum physics approach to 
conformal field theories, such categories appear. One can show that for rational 
conformal field theories on the circle, this category is in fact modular [M]- A 
modular tensor category is a special type of fusion categories that is highly non¬ 
degenerate. It is precisely such anyon models that are of interest to topological 
quantum computing. 

Many of the models that have such anyonic excitations are what is called 
topologically ordered. Topological order is a new type of order that does not fall 
into the Landau theory of spontaneous symmetry breaking, and there is no local 
order parameter distinguishing different phases. The examples known are not 
relativistic theories, and most of them are not described by a conformal field 
theory (at least, not directly). For example, the model that we will study here, 
Kitaev’s quantum double model, is a quantum spin model, defined on a lattice. 
Nevertheless, one could try to apply the ideas of Doplicher, Haag and Roberts 
to these systems, to study the superselection sectors and the properties of the 
anyons. This is indeed possible [MlH5] . 

This approach gives a strong mathematical founding to the study of anyons 
in topologically ordered systems. Not only does it allows us to borrow tech¬ 
niques from local quantum physics, it also opens up the way to the use of, for 
example, operator algebraic methods (see [48] for such an application). A rigor¬ 
ous mathematical framework may also be of use in this fast-moving field, where 
no consensus on the right definitions for fundamental concepts is reached. A 
precise setting where operator algebraic methods can be used may be the right 
setting to study important questions such as about the stability of such sys¬ 
tems: suppose that the Hamiltonian of the system is perturbed by a suitable 
perturbation, how much of the structure remains? One expects that the topo¬ 
logical nature of the system will preserve the interesting properties (as long as, 
for example, the Hamiltonian remains gapped), but this is something that one 
wants to prove rigorously (however, see e.g. [HKIl] for results in this direction). 

The remainder part of this contribution is outlined as follows. First, we 
discuss the basic idea behind topological order, and introduce the quantum 
double model in the setting of local quantum physics. Section [3] discusses the 
technical property of Haag duality in this context. We then come to the main 
part: an overview of the DHR-type analysis of the superselection sectors in 
the quantum double model. There we also relate it to the theory of (modular) 
tensor categories. This construction only works for the abelian quantum double 
models. Therefore, in the following section (which contains new results) we 
outline how one could use so-called amplimorphisms to study the non-abelian 
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case. Finally, in the last section we briefly comment on applications to quantum 
computing. 


2 Topological order 

Around the late ’80s, it was found that there are states of matter that do not 
fall into the Landau theory of spontaneous symmetry breaking. One of the 
first papers discussing this is [32]. Such states were called topologically ordered^ 
because some of their properties depend on the topology of the manifold on 
which the system is defined [62] . As a particular example one can think of two 
dimensional systems defined on a surface. For topologically ordered systems, 
the ground state degeneracy typically depends on the genus of the surface. 

Around the change of the millennium, quantum information theorists started 
to take an increased interest in topologically ordered systems. One of the main 
reasons is that it was realised that the topological properties of such systems 
might be useful for quantum computation. This was first suggested by Ki- 
taev m, who introduced a simple quantum system, the toric code, that could 
act as a quantum memory, storing a pair of qnbits. A qubit is the quantum 
analogue of a bit in a classical computer. Mathematically it is nothing but a 
copy of the Hilbert space C^, whose basis is usually denoted by |0) and |1). 
Although the basics of quantum computing are not difficult to explain, we refer 
to the standard textbook m for more information. In any case, the idea is 
to embed the two qubit Hilbert space into the physical Hilbert space 

H. We can then initialize our system in one of the states. The point of a good 
quantum memory is then that after some time, we should be able to recover the 
state (at least through the gathering of measurement statistics). In practice, 
however, quantum systems are not completely isolated, and influences from the 
environment might drive the system into a different state. 

This is where topologically ordered systems come into play. A nice feature of 
topologically ordered systems is that there is no loca0 order parameter. Hence, 
local observables cannot distinguish between different ground states, and on the 
other hand local operations cannot put the system from one ground state into 
another. Hence the idea is to encode information in the ground space, since 
this will be robust, at least against local perturbations. Of course, nature does 
not make it easy for us, and there are some more subtle points to being a good 
quantum memory. For example, because the information is stored non-locally, 
accessing it or acting on it also necessarily requires non-local operations, which 
may be difficult to implement. In addition, although local perturbations of 
the system do not destroy the information, this requires the use of an error 
correcting protocol. One should be able to do this fast enough, to prevent 
the errors from spreading out over the system to a non-local error, which does 
corrupt the information in the memory. At least for a wide class of 2D systems, 
it turns out that these systems by themselves are not good memories, so one 
indeed has to do some error correction mnEsunsj- 

Another interesting aspect of topologically ordered systems is that they gen¬ 
erally have anyonic excitations. That is, excitations of such systems behave 
like quasi-particles, with anyonic or braided statistics: the global state of the 
system changes non-trivially (that is, differently from just a sign change) under 

^For finite systems, “local” means small compared to the system size. 
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the interchange of two such quasi-particles. It is this property that we will focus 
on here. In particular, we will outline how the Doplicher-Haag-Roberts analysis 
of superselection sectors in algebraic quantum field theory can be translated 
to the setting of topologically ordered systems, and how this can be used to 
recover the statistics of the anyons. To illustrate how this works we consider 
Kitaev’s toric code, and more generally his quantum double model |37j . which 
is the prototypical example of a topologically ordered system. 

Since the toric code is relatively simple and nevertheless has many interesting 
features, one can use it as a testbed to see if ideas from algebraic quantum 
field theory can be applied to this class of models. This indeed turns out to 
work very well, and allows one to explicitly realize many of the fundamental 
concepts in the theory of superselection sectors. This shows that the concept 
can be transferred from relativistic theories to non-relativistic condensed matter 
systems. Although the explicit constructions we present here depend on the 
specific knowledge of the toric code, these features are common to a range of 
topologically ordered systems, and the theory can in principle be applied to 
them as well. 

To make the connection between the toric code and the theory of superselec¬ 
tion sectors, it is convenient to depart from the usual setting of Kitaev’s model 
and define it on an infinite 2D lattice, instead of on some compact surface of 
genus g. This is most conveniently done in the operator algebraic framework, 
were one assigns algebras of local observables to finite regions of space liiin]. 
In particular, there is a clear distinction between local and global observables, 
just like in local quantum physics. That is, one does not have to keep track of 
the system size. Concretely, we consider the lattice I?, and write T for the set of 
edges or bonds between them. We give each edge an orientation. For simplicity, 
all vertical edges are assumed to point upwards, the horizontal edges will point 
to the right (see figure [1]). Now let G be a finite dimensional group. Then to 
each edge we assign a “G-spin”. This means that there is a quantum system at 
each edge, with corresponding Hilbert space He '■= C[G], where the right hand 
side is the group algebra of G, seen as a Hilbert space in the natural way. 

The local operators that act on an edge e are *B(He) = M|g|(C). If A S 
'Pf(T), the set of finite subsets of T, then the local observables associated to A 
are defined as 21(A) := ^2 there is a natural inclusion of the 

corresponding algebras, by tensoring with the identity operator at the sites of A 2 
that are not in Ai. Hence we are in the familiar setting of local quantum physics, 
where we have a net A 21(A) of observables associated to bounded regions 
(cf. [24]). It is not difficult to see that this net is local, that is, [2l(Ai), 2l(A2)] = 

{0} if Ai n A 2 = 0. 

We then proceed in the same way as for relativistic systems: the local op¬ 
erators are defined as 2lioc = Uagp (r) 21(A). Note that 2t;oc is a =i=-algebra in a 
natural way, and that the standard operator norm on matrix algebras induces 
a G*-norm on 2l;oc- The quasilocal algebra 21 is the completion of 2lioc with 
respect to this norm. If A C T is an arbitrary subset (not necessarily finite), we 
then define 

—n- 

2t(A) = IJ 2t(A;) 

A/GP/CrnA) 

as the algebra of all observables that can be measured inside the region A. 

Finally, note that there is a natural translation symmetry on the system, 
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although we do not need this for our purposes, except for picking out a trans- 
lationally invariant ground state. This touches upon one of the fundamental 
differences between relativistic quantum field theory and discrete systems. For 
the discrete systems there is no Lorentz group or Poincare group, which play 
an important role in relativistic theories. Nevertheless, one can mimic some 
of these concepts using the translation symmetries mentioned, together with 
locality estimates for local observables evolved under the time evolution of the 
system [JS]. For our purposes these aspects play no role. 

2.1 The quantum double model 

So far the description has been completely general. To consider a specific system 
one has to specify the dynamics of the model. Note that because the model is 
infinite, the corresponding Hamiltonian generally is unbounded, and hence not 
an element of 21. Fortunately, in practice there is a lot more structure, and the 
dynamics are local in a suitable sense. More concretely, for each A £ T/(r) 
one can define a self-adjoint describing the interactions within that region. 
Heuristically, the dynamic evolution of an observable A is then obtained as 
at{A) := limA-j-oo where A —>• oo means that we take an increas¬ 

ing sequence of finite sets A that exhaust F. If the strength of the interaction 
decays fast enough, this expression converges and one obtains a strongly con¬ 
tinuous one-parameter group 1 1 —^ at of automorphisms m- 

Once the dynamics are defined one can talk about ground state. The most 
convenient way to do this is in terms of (generally unbounded) ^-derivations 
S, which are obtained as the generator of an automorphism group t at- 
In the cases of interest to us these are simply obtained as (the closure of) 
^ • ^loc t 21/oc; with 

(5(A) := i lim [i?A, A]. 

A—>-oo 

Because the interactions are local and of finite range in the cases of interest, this 
converges, and at (A) = e‘'^(A). A state wq on 21 is then called a ground state if 
—iuJo{A*d{A)) > 0 for all A £ D{S). This (at first sight) perhaps strange looking 
condition can be interpreted as a positivity of the energy condition. Such a state 
is automatically invariant, ujq o at = ujq-, so that in the GNS representation the 
dynamics are implemented by a strongly continuous group of unitaries 1U(t). 
By Stone’s theorem one obtains a Hamiltonian H, which can be shown to be 
positive using the condition mentioned above m- 

There is some extra notation that has to be introduced to define the dynam¬ 
ics in Kitaev’s quantum double model. A combination s = (v, f) of a vertex v 
and a choice of an adjacent face / is called a site. To each site we associate the 
following operators. Let g G G. Note that for a vertex v, there are four edges 
that start or end in u. A basis for the Hilbert space of this four sites is given 
by Isi) 0 • • • <8) |g 4 ), with gi £ G. The operator Af acts on this basis vector by 
multiplying gt from the left with g if the corresponding edge points away from 
the vertex v, and gt gets sent to gig~^ if the edge points inwards. If h £ G, a 
projection is defined as follows. First, list the edges around the face / in 
counter-clockwise order, starting in v. On the basis labelled by group elements 
of the Hilbert space corresponding to these four edges, B’f acts as the identity 
if a{gi) ■ • ■ <j{g 4 ) = h, and zero otherwise. Here a{g) = g ii the corresponding 
edge is in the same direction as the counter-clockwise labelling, and a{g) = g~^ 
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otherwise. Pictorially the operators can be depicted as follows: 


91 

94 

93 

919 ^ 

994 

993 


92 


929^^ 



93 


93 

Bf 

94 

91 

92 ^34 ^ ' 

94 

91 


Sometimes one says that projects on the states with flux h through the face 
/. Indeed, the structure is very similar to that of lattice gauge theory [53], and 
the quantum double model can in fact also be interpreted in such terms [37] . 

The operators Ag and Bg are also called star and plaquette operators respec¬ 
tively, and satisfy the following algebraic relations, which can easily be verified: 

AIA': = Af, BiB': = 5,,hBt A^B^i = Bf^~^A<i, (Af)* = Af\ (2.1) 

Star and plaquette operators acting on different sites s and s' always commute. 
At this point the name “quantum double” can be explained: the algebraic 
relations above are exactly those of the quantum double I?(G) of the Hopf 
algebra C[G] (see e.g. [33] for an introduction). That is, at each site there is 
an action of the Hopf algebra which acts via the star and plaquette operators. 
The representation theory of BiG) plays an important role in the analysis of 
the superselection sectors, as we will see below. 

We are now in a position to define the dynamics of the system. First, write 
Ag = Y^gfzc the average of the star operators, and Bg := B^. Then 

[As, Bg] = 0. The local Hamiltonians are then defined as 

i7A = -^As-^Hs, (2.2) 

sGA sGA 

where in the first sum the summation is over all sites s = (f,/), such that 
the edges starting or ending at v are contained in A. Similarly, the second 
summation is over all sites such that the edges of the face / are all contained in 
A. The case G = Z 2 is called the toric code. These local Hamiltonians generate 
a one-parameter group of automorphisms at as described above. It turns out 
that there is a unique translationally invariant ground state [I] |46] H 

Proposition 2.1. There is a unique translationally invariant ground state 
loq. This state is pure and is the unique state on 21 which satisfies a;o(As) = 
ujQ{Bg) = I for all star and plaquette operators Ag and Bg. 

This ground state will be the starting point from the analysis: the differ¬ 
ent superselection sectors will be realized, roughly speaking, by creating single 

^Bruno Nachtergaele pointed out that the remark in 1451 is in fact false, and there are ad¬ 
ditional (non-translationally invariant) ground states. Indeed, the charged states constructed 
in m have dynamics implemented by a positive Hamiltonian. 


6 









Figure 1: The lattice describing the system, together with a ribbon between two 
sites. Note that the ribbon has an orientation, pointing from start to end. 


excitations of this ground state. The GNS representation corresponding to the 
ground state representation from the proposition will be denoted by (tto, fl, 'Hq). 
Since 21 is a UHF (and hence simple) algebra, ttq is automatically injective, and 
we will often identify 7ro(A) with A. 

To understand how this works it is necessary to understand what we mean 
by an excitation, and how these can be obtained. From the proposition above 
it follows that in the ground state representation, = n. We can 

interpret this as constraints, and violating some of these constraints carry an 
energy penalty, according to the local Hamiltonians (12.21) . We will interpret the 
violation of such a constraint as a (quasi-)particle sitting at the site s (if the 
state is an eigenstate of i?A)lll So the excitations live at sites of the lattice. 

The excitations can be obtained by acting with so-called ribbon operators 
on the ground state. We only list their main properties here. A more complete 
treatment and proofs can be found in A ribbon is, roughly speaking, a 
continuous path of triangles (see figure [T] for an example). To each such ribbon 
^ and pair g, h of group elements one can assign a ribbon operator . These 
operators act only on the edges that are part of the ribbon (or cross any of the 
triangles). They satisfy the following algebraic relations: 

TT'^ijpi T^h2,92 r Tphih 2 ,gi { _ ph ^ ,g \ ^ _ T 

~ ^ 91 , 92 -^^ 5 , 2 ^^^ — -'• 

g^G 

Note that the ribbons carry an orientation. This allows us to talk about the 
starting and ending sites of the ribbon. Let us denote them by si and S 2 for the 
moment. An important property of the ribbon operators is that they commute 


^Note that we disregard any momentum variables. The pairs of excitations form bound 
states, so the excitations mentioned here are not quite the single-particle excitations one 
encounters in scattering theory. 
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with any star or plaquette operator, except for those at the ending sites of the 
ribbon. There we have the following commutation relations: 


_ j^khk ^ ^kg j^k j^k ph,g j^h,g j^kh 


Similarly, for the ending site 


Ak Tph^g j^h.gk ^ jfc rik rph^g j^h,g j^g ^gk 

^S2> ^S2 


(2.3) 


(2.4) 


From these commutation relations it follows that the vector state violates 

some of the constraints the Hamiltonian gives for the ground state, namely 
precisely those at the start- and ending sites. Hence this vector state can be 
thought of as an excited state of the system. These excited states satisfy the 
following important property, whose proof we will omit (see e.g. [3 [231 for the 
proof). 

Lemma 2.2. Let and ^2 be two ribbons with the same starting and ending 
sites. Then 

In other words, the state does not depend on the ribbon itself, only on the 
endpoints. This already shows the topological nature of the system. 

Now suppose that ^ is a non-trivial ribbon. It can then be written as the 
concatenation of ribbons and ^ 2 - The ribbon operator ^ is related to the 
smaller ribbons by means of the following recursion relation (which indeed can 
also be used to recursively define the ribbon operators): 


Tph,g _ \ Tph,k T^khk,kg 

k&G 


(2.5) 


This feature will be important later. 

The basis F^’^ of ribbon operators is not always the most convenient one. 
Recall that for each site s there is an action of T>{G) on the Hilbert space. One 
can now consider the vector space V = spa,n{F^’^Ll : h,g G G} for a hxed 
ribbon By equations (ESI) this space has a natural interpretation as a left 
T’(G)-module, because = fl and = Sh,e^. Similarly, equations (12.411 
also give it the structure of a I>(G)-module, induced by the action of the star 
and plaquette operators at the ending site0 This suggests that it may be useful 
to decompose this space into irreducible subspaces, and find a new basis of the 
ribbon operators accordingly. This indeed turns out to be a good idea. 

The representation theory of T>(G) is well understood: constructions of all 
irreducible representations can be found in B mi- They are in one-to-one 
correspondence with pairs {C,p), where G is a conjugacy class of G and p an 
irreducible representation of Zg{c), where Zg{c) is the centralizer in G of a 
representative c of the conjugacy class. It turns out to be convenient to label 
the elements of the conjugacy class in a particular way (compare [5])- First of 
all, let G be a conjugacy class of G. Choose a representative r G G, and let 
Zair) be the centraliser of r in G. We label the elements of G by ci,... ,c„, 
where n = |G|. Then there are Qi such that Ci = Qirq^., where we used the 
notation g for the inverse of g, to improve readability. The set {qi} is denoted 


^This is in fact the contragradient module of the representation at the starting site. 



by Qc- Note that each g € G can be uniquely written as g = giU for some 
Qi e Qc and n e Zc{r). 

With this notation we can choose a new basis of the ribbon operators asso¬ 
ciated to a ribbon Suppose that p is a unitary representation of Zcir). We 
regard each p[g) as a unitary matrix, and write p{g)jj' for the corresponding ma¬ 
trix elements in the standard basis. Let iQ' = 1,... n and j, j' = 1,... dim(p). 
We then define 


,j,j' _ 


H Pjj'i9)F^ 




geZcir) 


As C runs over all conjugacy classes of G, and p runs over the corresponding 
irreducible representations of the centralisers, these operators form a basis of the 
space spanned by . We refer to [5] for a proof. In essence, the point is that 
the space of operators is decomposed into subspaces transforming according to 
some irreducible representation of LI(G). For convenience we sometimes drop 
the notation for C and p when these are implied by the context, and write 
I = (i,j), J = (*',/) for the pairs of indices: . 

If the group G is abelian the notation can be simplified, since each conjugacy 
class consists of exactly one element, and clearly the centralizer Zc{r) is always 
equal to G. Hence we can label the basis by a pair (w, c), where w is a character 
of G and c € G. It can be checked that in that case the corresponding ribbon 
operators F^’^ are unitaries and that = p^x,cd^ There are many more 

useful relations, describing for example the commutation relations between two 
ribbon operators acting on crossing ribbons, these can be found in 


3 Haag duality 

Just as in the Doplicher-Haag-Roberts theory of superselection sectors, Haag 
duality plays a very useful role. Let us first recall what Haag duality actually is: 
it is an commutation property of von Neumann algebras generated by algebras 
of observables in the ground state representation. More precisely, let A be a 
cone-like region^ Then 7r(2t(A))" is the von Neumann algebra generated by all 
(quasi-)local observables localised in A. Note that, by locality, these observables 
commute with any local observable localised outside the cone, that is, in A*^. 
In other words, 7ro(2t(A))" C 7ro(2t(A‘^)', where A° is the complement of A in T. 
Haag duality says that these sets are actually equal. It means that we cannot 
add operators to the cone algebra tto( 21(A))" without violating locality. This 
property is fulfilled for Kitaev’s model for abelian groups G: 

Theorem 3.1. Let G be a finite abelian group and let K be a cone. Write ttq 
for the translational invariant ground state representation of the quantum double 
model for G. Then Haag duality holds: 

TTO (21(A))" = 7ro(2t(A^))', 

where the prime denotes the commutant in the set of all hounded operators. 


^The precise shape is not so important, but see m for a precise definition. Heuristically, 
one can take a point in the lattice, draw two semi-infinite lines from this point, and consider 
all edges that either intersect these lines, or lie in the convex set bounded by the lines. 
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This theorem was first proven in m for the case of G = Z 2 (the toric code) 
and later extended to all finite abelian groups in ^5] . 

The proof of the theorem depends on a good knowledge of the pairs of 
excitations and the operators that generate them. This makes it possible to 
find a convenient description of the Hilbert space of the system, and the Hilbert 
space 'Ha describing pairs of excitations localised in a fixed cone A. With this 
description it is possible to demonstrate that the self-adjoint parts of certain 
algebras (the restriction to 'Ha of the algebras of observables localised in A, and 
the algebra of observables localised in A'^), when acting on the cyclic GNS vector 
generate a dense subset of 'Ha- It is known that this is related to commutation 
properties of algebras m. and this allows one to conclude Haag duality. 

Much of the proof boils down to a good understanding of the representa¬ 
tion theory of T’(G), the quantum double of G, since the operators that create 
pairs of excitations correspond to irreducible representations of the Hopf alge¬ 
bra. Their behaviour under the interchange of two such operators (or on the 
representation theory side, the braiding of the category of representations) also 
plays a role. These properties are well-known and studied for a wider class 
of quantum doubles, particularly those associated to Hopf-* algebras that are 
quasi-triangular (which means a braiding can be defined). As Kitaev already 
remarked in his paper the quantum double model can also be defined for 
such Hopf algebras. The following conjecture therefore seems very natural: 

Conjecture 3.2. Consider a quasi-triangular Hopf-* algebra H. Then the 
corresponding Kitaev model on the plane satisfies Haag duality for cones. 

The main difficulty in proving this is that the combinatorics get much more 
involved. In particular, in the case of abelian algebras all irreducible repre¬ 
sentations are one dimensional. For non-abelian algebras one has to consider 
multiplets of operators that create pairs of excitations, transforming according 
to some irreducible representation of the Hopf algebra V{H). We will see some 
of the consequences later on, when we discuss non-abelian theories. 

It should be noted that many of the constructions we present below do not 
depend on Haag duality. The only place where it is used is in going from a 
representation that satisfies the selection criterion to a localised endomorphism. 
Representatives of these endomorphisms (and intertwiners between them that 
change the localisation region) can be constructed without an appeal to Haag 
duality. However, it is then not possible to conclude that each representative of 
the equivalence class is indeed given by a localised endomorphism. 

4 Superselection theory for abelian models 

It is well-known that the superposition principle of quantum mechanics does 
not hold unrestrictedly. A familiar example is that of bosons and fermions m- 
Consider the vector if = ■^('*A/ + f’b), where ■i/'/ is a single-particle fermionic 
state, while tpf is bosonic. Under a rotation of the system by 360 degrees, the 
fermionic state will acquire a minus sign, while the bosonic part is unchanged. 
Physically, however, the states are indistinguishable, and will lead to the same 
expectation values. In general, one can argue that that the physical Hilbert 
space can be decomposed into different sectors corresponding to the different 
types of charges in the system. This is called a superselection rule. 
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On a mathematical level, superselection sectors arise because there are in¬ 
equivalent representations of the observable algebra 21. A (7*-algebra has many 
inequivalent representations in general, hence one should somehow select the 
physically relevant representations. There are several conditions one might im¬ 
pose. For example, in relativistic theories a natural condition is to look at 
representations that are covariant with respect translations, such that the spec¬ 
trum of the generators (that is, the momentum) is contained in the forward light 
cone. Doplicher, Haag and Roberts proposed to look at those representations 
that, roughly speaking, look like the vacuum representation in the spacelike 
complement of a double cone, the intersection of a forward and backward light- 
cone m- Although this certainly does not cover all physical systems of interests 
(it excludes, for example, electromagnetic charges [H]), it is nevertheless a use¬ 
ful criterion, and can for example also be applied to conformal theories on the 
circle P^] . 

What Doplicher, Haag and Roberts (DHR) showed [Tallin] (see also 
is that these representations can be studied in a systematic way. In this way, one 
can learn something about the statistics of the different charges (i.e., if they are 
bosons or fermions in space-times of dimensions > 2-1-1), or how they behave 
under composition of two charges (so-called fusion rules). Mathematically, it 
amounts to studying the structure of the equivalence classes satisfying the se¬ 
lection criterion as a braided fusion category^ which is a type of tensor category. 
Some of the essential steps will be outlined below, but a more thorough intro¬ 
duction can be found in, e.g., [501 [SI]. The study of tensor categories is a whole 
field on its own (an introduction can be found in [42] 1. but familiarity with the 
main terminology is not strictly necessary for our purposes here. 

In the remainder of this section G will be a finite abelian group. We will 
see how one can construct different superselection sectors for Kitaev’s quantum 
double model, and do a Doplicher-Haag-Roberts type of analysis to recover the 
properties of the different charges. The results outlined here have been obtained 
in [231145] ■ to which the reader is referred for details. It is perhaps somewhat 
surprising to see that this theory, which was originally developed for relativistic 
quantum systems, can be applied so successfully to discrete lattice quantum 
spin systems, even though there are fundamental technical differences. 

4.1 Localized representations 

As mentioned, the different superselection sectors are identified as equivalence 
classes of representations of the observable algebra, satisfying additional selec¬ 
tion criteria. Here we take the ground state representation tto, corresponding 
to the translationally invariant ground state, as a reference representation, and 
look at all representations that look like ttq when considering observables out¬ 
side a cone-like region A. The reason for this will become clear later, but at 
this point we mention that this is similar to the work of Buchholz and Freden- 
hagen, who show that in relativistic theories, massive particles can be localised 
in spacelike cones usmn, leading to a similar criterion. We will write C for the 
set of cones. 

Definition 4.1. A representation tt o/2l is called localizable if for each A G C 
we have 

TTO r2l(A^)^7r r2l(A=), (4.1) 


II 


where = denotes unitary equivalence of representations and f means that we 
restrict the representation to the subalgebra 2t(A'^). 

An equivalence class of representations satisfying the criterion is called a 
(superselection) sector or simply a charge. 

To proceed in the analysis of these representations, we have to pass from 
representations to endomorphisms of the observable algebra. This can be done 
with the help of Haag duality. Indeed, let A S £. Then from equation (ED 
there is a unitary Va setting up the equivalence with ttq for the subalgebra 
2t(A'^). Define p(A) = Va_'k{A)V^ and let A 2 G £ contain A. Let A G 2l(A2) and 
B G 2t(A2). Note that p(B) = 7ro(B). Moreover, by locality, 

p(AB) = p{A)tto{B) = p{BA) = ■Ko{B)p{A), 

hence by Haag duality p{A) G 7ro(2l(A2))' = 7ro(2l(A2))". Since ttq is faithful, 
we can identify 21 with 7ro(2l) and regard p as an endomorphism of 2lH The 
endomorphism is called localised in A. It is also transportable: if A' G £, 
there is a unitary V and endomorphism p' (localised in A') such that Vp{A) = 
p'{A)V for all A G 21. Such an operator V (not necessarily unitary) is called an 
intertwiner from p to p'. If H is unitary we will also call it a charge transporter, 
since it moves a charge from one cone to another. Using Haag duality one can 
show that in fact V G 21(A) for suitable A € C (where A should contain the 
localization regions of both p and p'). There is in fact a 1-1 correspondence 
between superselection sectors and equivalence classes of these localised and 
transportable endomorphisms. From now on we will work with the latter. 

4.2 Localized sectors in the quantum double model 

Recall that the ribbon operators create a pair of excitations (or charges). 
We are however interested in the properties of a single charge. Since we are in an 
infinite system, it is possible to create a pair of excitations, and move one of them 
to infinity. As may already be anticipated from the discussion in Section 12.11 
the different charges are in 1-1 correspondence with pairs (w, c) of a character w 
and group element c. The ribbon operator then creates the corresponding 
charge at the beginning of and a conjugate charge at the other end. On the 
level of the observables, this means that we map A Hence 

to implement the idea of moving one charge to infinity, one can choose a semi¬ 
infinite ribbon which for simplicity we assume to be completely inside some 
A G £, and write for the (finite!) ribbon consisting of the first n parts. Then 
we define 

a{A) = lim F;"’"A(F;"’")*. (4.2) 

Considering the dense set of local observables, and using the decomposition 
rule (1^ . it is not so difficult to show that this expression converges, and 
defines an automorphism of 21. It describes how the observables change in the 
presence of a single charge (a;, c) in the background. 

^Strictly speaking, this is only true if 2l(A2)^^ C 21, which in general is not the case for 
unbounded regions A 2 . One can however solve this by passing to a larger algebra 21"^“ D 21 
and extend p to a proper endomorphism of 21-^“ |15l Sect. 4]. 

^This actually corresponds to creating charges with , which turns out to be slightly 

more convenient. 
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The map a defined above leads us to an example of a representation that 
satisfies the selection criterion, by defining tt = ttq o a. Clearly, if A is any 
cone containing the ribbon, a{A) = A for all A G 2t(A'^) by locality, so a is 
localised. It is also transportable. There is a neat way of seeing this: because 
the vector state depends only on the endpoints of the path, it follows 

that two for automorphisms ai and defined in terms of the same charge, 
but different ribbons (with the same fixed endpoints), the states loq o cei and 
Wo o 02 coincide. In addition, note that both triples (tto o ai,r2,Ho) are GNS 
triples for this state, so that by the uniqueness of the GNS construction, ttq o ai 
must be unitarily equivalent to ttq o 0 : 2 . The argument can be extended if 
the endpoints of the ribbon do not coincide, by conjugating with a suitable 
ribbon operator. It follows that the automorphisms are transportable. The 
corresponding intertwiners can be shown to be in 7ro(2t(A))" for a suitable cone 
A using Haag duality, but for the quantum double model an explicit construction 
of a net converging in the weak operator topology to an intertwiner is also 
possible. 

The construction above gives for each pair (a;,c) an equivalence class of 
localised and transportable automorphisms, but it is not clear yet that these 
classes are distinct. This does turn out to be the case. The key idea in showing 
this is by considering an analogue of Wilson loops, which measure the charge in 
the region enclosed by the loop. Such operators can be obtained by considering 
ribbons with the same start and ending sites. In this way one can construct a 
charge measurement operator in an arbitrarily large region, that have expecta¬ 
tion value I in the state wq o a, where a has charge (uj, c), and zero expectation 
value in states obtained from a different pair It follows that the cor¬ 

responding representations ttq o a must be inequivalent. The discussion can be 
summarized in the following theorem: 

Theorem 4.2. Let G be an abelian group. For each pair (w, c) of a char¬ 
acter of G and an element c G G, there is an equivalence class of localised 
and transportable automorphisms. If ai and 02 are such automorphisms, then 
ttq o ai = TTo o 02 if and only if the belong to the same charge class (w, c). 

This result can also be understood as an instance of charge conservation. 
The total charge of the state can be obtained by finding all excitations w (resp. 
c) in the state, and multiply all of them together. Note that this is possible 
because the dual of an abelian group is a group again. The ribbon operators 
span a dense subset of the quasilocal algebra 21. Because the two distinct charges 
at the end of the ribbons transform according to conjugate representations, it 
follows that one cannot change the total charge of the system just by local 
operations. In particular, one cannot go from a charged sector with a certain 
total charge to another sector with a different total charge by acting with local 
operators, so that indeed we have inequivalent representations of 21. 

4.3 Braiding and fusion 

So far we have constructed endomorphisms that describe single charges localised 
in cones, and charge transporters or intertwiners that can move the charges 
around. There is however much more structure, and this is where it becomes 
essential that we have endomorphisms rather than representations: in contrast 
to representations, endomorphisms can be composed. If pi and p 2 are two 
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localised endomorphisms, we define pi ® P 2 {A) = pi o p 2 {A). The interpretation 
is that we first create a charge p 2 , and then add a charge pi to the system. This 
operation is called fusion. Note that pi o p 2 is localised again, more particularly 
it is localised in any cone that contains the localization regions of pi and p 2 . It 
is also transportable again. This can be seen by defining a product operation 
for intertwiners as well: if Spi{A) = p 2 {A)S and Tai{A) = a 2 {A)T for all 
A € 21, we define S ® T := Spi{T). Then an easy calculation shows that 
S ® T intertwines pi ® cti and p 2 ® (72 ■ There is however a slight technical 
issue. The intertwiners S are in general only elements of weak closures of the 
form TTo(21(A))" for some cone A (this follows from Haag duality), while the 
endomorphisms are a priori only defined on 21. It is however possible to extend 
the endomorphisms to a slightly larger algebra 2l^“, where A^ is some fixed 
auxiliary cone, that does contain the intertwiners. Since this is a minor technical 
point, we ignore the issue here, and refer to [13 sa for technical details. 

Even if pi and p 2 are irreducible, their composition need not be irreducible 
any more. A natural question therefore is if it is possible to decompose the 
composition into irreducibles again. The fusion rules give this decomposition. 
That is, if pi and pj are irreducible endomorphisms, there are integers A^*-, where 
k runs over a set of representatives of all irreducible localised and transportable 
endomorphisms, such that 



(4.3) 


k 


The sum operation is a direct sum, which can be described in terms of intertwin¬ 
ers. This is very similar to representation theory of finite (or compact) groups: 
the tensor product of two irreducible representations can be decomposed again 
as the direct sum of finitely many irreducible representations 0 Note that the 
identical endomorphism t(A) = A acts as a unit for the fusion operation. Fi¬ 
nally, a charge p has a dual or conjugate p if there is such a localised and 
transportable p with = 1. The existence of conjugates is not automatic, 
and has to be either proven or taken as an assumption. Physically, it is related 
to the existence of (charge) anti-particles. 

In algebraic quantum field theory, it is usually shown that such a decomposi¬ 
tion exists with the help of a technical property. Property B [S] , that essentially 
allows any projection in the local algebra to be written in the form P = WW* 
for some isometry W, localised in a slightly bigger region. In the present sit¬ 
uation of spin systems, one could go about by showing that direct sums can 
be constructed, and by explicitly verifying the fusion rules (14.31) [IS]. This is 
particularly easy in abelian models, where there always is a unique fusion out¬ 
come. That is, Njj is equal to one for exactly one value of k, and zero otherwise. 
Symbolically, for the abelian quantum double model the fusion rules are 


(a;i,ci) 0 (a; 2 ,C 2 ) = (wia; 2 , C 1 C 2 ). 


This can be seen by considering, without loss of generality, a single ribbon f, and 
define the automorphisms ai and 02 corresponding to the choices of charges. 


^Indeed, if the theory contains only bosonic sectors, the localised and transportable en- 
domorphisms are in one-to-one correspondence with the representations of a compact global 
symmetry group G mi- 
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Using the multiplication rule for two ribbon operators in abelian models, men¬ 
tioned at the end of Section [O the result follows. 

The final piece of structure that we discuss here is that of braiding. This 
is related to the statistics of identical particles, that is, their behaviour under 
interchange. This amounts to the study of the relation between pi ® p 2 and 
P 2 ® Pi for general localised endomorphisms pi and p 2 - An intertwiner £pi,p 2 
relating the two can be constructed explicitly by, quite literally, moving p 2 
around (or, pi of course, in similar way). More precisely, one can choose a 
localization region A that is disjoint from the localization regions of both pi and 
P 2 . By transportability there is a unitary U and a endomorphism p localised in 
A such that U P 2 {A) = p{A)U. Because pi®p = p®pi since their localization 
regions are disjoint, it follows that epi,p 2 := {U* ® I){I ® U) = U*pi{U) is an 
intertwiner from pi o p 2 to p 2 o pi. 

In three or more spatial dimensions, the definition of epi,p 2 is completely 
independent of the choices made, and one can show that £pi,p 2 £p 2 ,pi = Note 
that this is the situation of ordinary bosons and fermions: moving one particle 
around the other doesn’t change the system. One can show that the operators 
Ep^p indeed induce a representation of the permutation group, which interchanges 
charged excitations |19] . In lower dimensional space times, things get more 
interesting. It is no longer true that there is a unique choice of £pi,p 2 [ISl [551 
128] . The reason is that one cannot continuously move the cone A around, but 
instead has to choose between either “left” or “right”. This can be defined 
unambiguously by choosing an auxiliary cone A^ as above. To get a consistent 
definition for the operators £pi,p 2 (for example, to make sure that the different 
ways to go from {p^a^r to r^p^a coincide), one has to choose one of 
the two alternatives, and stick to that. A consequence of this ambiguity is that 
it is no longer true that £pi,p 2 £p 2 ,pi = I- In. addition, one does not obtain a 
representation of the permutation group any more, but rather one of the braid 
group. Consequently, such charges are said to have braid statistics. If this 
representation is abelian, one calls the charges (abelian) anyons (because they 
can pick up “any” phase under interchanged If the representation is non-abelian, 
one usually speaks of non-abelian anyonsjj 

The charges of the quantum double model are abelian anyons. This can be 
verified directly, by calculating the operators epi,p 2 s-nd £p 2 .piJ together with 
the explicit construction of charge intertwiners, mentioned above. In the end 
the calculation boils down to commutation relations of the ribbon operators: 

= a;(d)cr(c)F^’‘^FK’“ for two crossing ribbons and ^ 2 - In the end 
one finds that 

^(d)u'(c), 

showing that the charges are indeed abelian anyons. 

4.4 The category of localised endomorphisms 

The set of localised endomorphisms has a very rich structure, of which we have 
given some examples above. Mathematically, it has the structure of a braided 


®In the local quantum physics literature the name plektons is used. This name does not 
seem to have caught on outside of that community: the name non-abelian anyons prevails for 
example in the field of topological quantum computing, and this is why we adhere to that 
name here. 
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unitary fusion category^^ The theory of such categories is quite rich, and an 
active field of study, in particular considering the classification of such categories. 
There are by now a few texts that provide an accessible entry to the literature. 
We mention for example [ssiEiiiiiini], each of which has a different focus. 

We will denote this category by A. Its objects are the cone localised and 
transportable endomorphisms of 21. The morphisms are intertwiners: T S 
Hom(p, cr) if Tp{A) = < 7 {A)T for all A G ‘Ql. Note that if T € Hom(/9, cr), 
then T* G Hom((T, p), where T* is the adjoint of T. The isomorphisms in this 
category are unitary operators, so that the we indeed obtain equivalence classes 
of localised and transportable endomorphisms. An object is irreducible, by def¬ 
inition, if Hom(p, p) = C. Note that this is just another way of stating that the 
commutant of p(2l) is trivial, hence this coincides with the usual notion of an 
irreducible representation. 

In the previous section we defined a tensor product. This is actually a tensor 
product in the category, where the trivial endomorphism is the tensor unit. It is 
in fact a strict tensor category: associativity of the tensor product holds on the 
nose, not only up to isomorphism as is often the case in category theory. For 
each pair of objects pi and p2, the definition of the braid operator £pi,p2 gives 
an isomorphism £pi,p2 G Hom(pi, p 2 ). This assignment is in fact natural in both 
variables, and satisfies the so-called “braid equations”, which give consistency 
conditions. This makes A into a braided tensor category. It is symmetric if 
^Pi’.P 2 = ^P 2 .pi foi' all objects pi,p 2 - 

Direct sums can be described in the category theory language as well. If 
Pi and p 2 are objects in A, an object pi © p 2 is a direct sum if there are 
Vi G IIom(pi, pi®p2) with ^ ■ ViV* = I and V*Vj = In algebraic quantum 
field theory the existence of such direct sums follows from Property B discussed 
above. In the case of the quantum double model we can explicitly construct 
the direct sums. Let A be a cone and suppose that pi and p2 are localised in 
A. Adapting the arguments of [3S1 [31], it follows that T^a := tto( 21(A))" is an 
infinite factor, and hence there are isometries Vi G TZa such that ViV* = I 
and ViV* = Sijj. We can then define 

Pi © P2(A) := Vipi{A)Vf + V2P2{A)Vf. 

Because of locality it follows that pi ©p 2 is localised in A, and a straightforward 
calculation shows that this indeed is a direct sum in the categorical sense. Note 
that the direct sum is only defined up to unitary equivalence. This also explains 
how the summation in the fusion rules has to be understood, namely as a direct 
sum in this category. Note that in the abelian models we consider here, there 
is only one non-zero fusion coefficient, so the direct sum construction is not 
necessary. 

Finally, there is the notion of duals. Physically, these can be interpreted as 
anti-charges. If p is an irreducible object of A, a dual (or conjugate) is a p S A 
such that the trivial endomorphism i of 21 appears exactly once in the direct 
sum composition of p(8>p and of p(8>p. Note that this coincides with the notion 
of a conjugate in terms of the fusion coefficients mentioned before. This 
implies that there is an isometry R G Hom(i, p © p). We assume all objects 

Strictly speaking, “fusion” implies that there are only finitely many equivalence classes 
of irreducible objects. This is the case for the quantum double models, but need not be true 
in general. Dropping that condition does not make a difference in many cases. 
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have conjugates. In this case, it follows that all Hom-sets are finite-dimensional 
vector spaces (over C). It follows from the fusion rules that for the abelian 
quantum double model, the conjugate of (w,c) is (oJ, c“^). 

So far we have constructed different superselection sectors of the abelian 
quantum double model, and studied their properties. The question remains if 
there are perhaps additional sectors, that we have so far overlooked. In other 
words, are there perhaps representations that satisfy Definition 14.11 but are 
not equivalent to one of the representations constructed so far? This question 
can be answered by adapting techniques from rational conformal nets |d4) . In 
particular, one can consider two disjoint cones Ai and A 2 , and look at the 
von Neumann algebra 7ro(2l(Ai U A 2 ))". While Haag duality holds for a single 
cone, for a pair of cones it generally does not hold anymore. This means that 
7ro(2l(Ai U A 2 ))" C 7ro(2t((Ai U K2Y)" in general is a proper inclusion (it is 
an irreducible subfactor, in fact). Essentially, the bigger algebra also contains 
the intertwiners or charge transporters that move a charge from one cone to 
the other. Hence by studying how much bigger the algebra is, one can learn 
something about the number of sectors in the theory. This relative size can 
be quantified by the Jones-Kosaki-Longo index [7ro(2t((Ai U K 2 Y)" '■ 7ro(2l(Ai U 
A 2 ))"], which in turn gives an upper bound on the number of sectors |48] . For 
the toric code this number is four, and hence it follows that we have indeed 
construct all sectors: the ground state sector, and the ones corresponding to 
the pairs (t,^), (ct, e), (cr, g), where a is the sign representation of Z 2 and g is 
the only non-trivial element in the group. This leads to the following conclusion: 

Theorem 4.3. The category of localised and transportable endomorphisms of 
the toric code is equivalent as a braided fusion category to category RepyT>(G') 
of finite dimensional representations of the quantum double 2 ?(Z 2 ). 

It should be rather straightforward (but tedious) to extend this result to all 
abelian quantum double models. 

There is one aspect that we have not mentioned so far. The category of 
the quantum double model is modular. Modular tensor categories are fusion 
categories with an additional non-degeneracy property: it is called modular if it 
has trivial braided centre. That is, suppose that p is irreducible, and ep^cr£a,p = I 
for all irreducible a. If this implies that p is the trivial object, then the category 
is modular [41] [55]. In this sense, modular categories are “as far away” as 
possible from being symmetric, where £p^a£cr,p is always the identity. 

There is an equivalent condition that is related to Verlinde’s matrix S in 
conformal field theory m- In the categorical setting it can be defined as follows. 
One can show that the duality allows one to define a trace on morphisms in the 
category. A matrix S is than defined by having entries Sij = trYp.^p^Sp^^pY, 
where pi is a set of representatives of the irreducible objects. A category is 
modular if and only if this matrix S is invertible. It also allows us to explain 
the name modular: together with a matrix T that can also be canonically 
defined, the matrices S and T induce a projective representation of the modular 
group SL 2 (Z). In modular categories it turns out that there is in fact a relation 
between the matrix S and the fusion coefficients, given by the Verlinde rule [60j : 

Arfc _ SjjSjrSkr 

~ 2 -^ C2 ’ 
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where the label 1 stands for the trivial object. This shows that that the structure 
of a modular tensor category is quite rigid, and not any given fusion rule can 
be realized in some modular category. 

Modular tensor categories can be realized as topological quantum field the¬ 
ories (TQFTs), see for example [59]. One can therefore think, in a sense, of 
the type of topologically ordered systems that we have considered here as sys¬ 
tems that in the low energy limit can be described by a TQFT. It should be 
noted than in general these are really effective theories, giving an effective de¬ 
scription of oft-times very complex and poorly understood condensed matter 
systems. They also play a fundamental role in the field of topological quantum 
computing |61] . on which we comment briefly in the last section. 


5 Extension to non-abelian models 

So far we have only considered abelian models. A natural question is if the meth¬ 
ods can be extended to non-abelian quantum spin models. We again consider 
the quantum double model, but now for a non-abelian (but still finite) group 
G. In that case, there still is a unique translational invariant ground state luq. 
Many of the proofs, however, do not directly carry over from the abelian case. 
The reason for this will become clear below, but the underlying difficulty is that 
the irreducible representations of T){G) are no longer all one-dimensional. In 
particular, it is not clear how one could construct endomorphisms describing 
these “non-abelian” charges. 

A similar problem appeared in one-dimensional spin chains, with compactly 
localised charges. There the problem is that the algebra 21(A) of such a local¬ 
ization region is finite dimensional, namely the tensor product of finitely many 
matrix algebras. All endomorphisms of this algebras are in fact automorphisms, 
and one can show that these cannot have non-abelian statistics. This problem 
can be circumvented by the methods of Szlachanyi and Vecsernyes jSS], and 
Nill and Szlachanyi m- Instead of looking at endomorphisms, they look at 
amplimorphisms^ that is, morphisms y : 21 —>■ Mn(%) for some integer n. One 
can then do a study of these endomorphisms in the spirit of the DHR theory. 

In the present case we are interested in two dimensional systems with lo¬ 
calization in cones. The cone algebras are certainly not finite dimensional, so 
the obstruction of the only endomorphisms being automorphisms does not play 
a role herell^ Nevertheless, as mentioned it is not easy to explicitly construct 
examples of the endomorphisms. This is where the amplimorphisms come in. 
The strategy is to mimic the amplimorphism construction in [58] in the context 
of cone-localised charges, to construct representatives of the different charge 
classes, and will then show how we can go back to the usual setting of cone- 
localised endomorphisms. The amplimorphism description is much more ex¬ 
plicit, making it possible to explicitly calculate intertwiners, fusion rules, et 
cetera. Here we will mainly restrict to the construction of the different sectors. 
Combining the techniques developed in the abelian case with the amplimor¬ 
phism results in [^ should enable one to completely solve the model. 


be a bit more precise, one actually needs that ttq (21(A))" is not a factor of Type I, 
that is, not of the form 33 (H) for some Hilbert space H. For the quantum double models that 
is the case. 
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The ribbon operators again play a fundamental role. We will use the same 
notation as in Section mi e.g. r will always be a fixed representative in a 
conjugacy class. Let ^ be a fixed ribbon, and recall that one can choose a basis 
of the ribbon operators acting on ^ in terms of the irreducible representations of 
T>{G). For a fixed representation there is a multiplet of ribbon operators , 
where I = and J = (jijjb)- The first indices run over the elements of 

the corresponding conjugacy class C, while the second runs over the dimension 
of the group. These multiplets satisfy certain completeness relations, namely 

= W, (5.1) 

I J 

where the summations are over all pairs (ii, 12 ) and the I on the right hand side 
is the unit of 21. These equations can be verified by a simple calculation using 
the definitions and the algebraic relations for . In these F^'^ are called 
irreducible and complete multiplets, with the difference that here we have not 
defined an action 7 ^ of V{G) acting on these multiplets. Nevertheless, they do 
transform under an irreducible representation of 'D{G), cf. equation (B69) of [5]. 

Recall that for the abelian model, a ribbon operator acting on a large rib¬ 
bon is just the product of two ribbon operators acting on smaller ribbons, but 
related to the same irreducible representation. If the irreducible representation 
is not one-dimensional any more, this is no longer true. Nevertheless, a suitable 
analogue of equation (j2.5l) still holds: 

Lemma 5.1. Choose a pair {C,p) of a conjugacy class and an irreducible rep¬ 
resentation of the centralizer Zcir), where r is as explained in Section \2. 1\ Let 
f, = ^ 1^2 be a ribbon that is decomposed into two ribbons. The correspond¬ 
ing multiplets will be denoted by F^'’ and F^^ respectively. Then we have the 
following relation 

F" = T. Pc (5.2) 

K 

where the sum is over all pairs K = (fci, ^ 2 ) 

Proof. First write out the right hand side of equation (15.21) in terms of the 
elementary ribbon operators F^j^, where we set / = ( 11 , 12 ), and similarly for J 
and K: 


E 


|C| dim(p) 

_ / x_ ,,.Ci^,qi^gq^ Ck^,qk^hq 


g,h£Za{r) ki=l k2 = l 


Since p is a representation, the summation over ^2 yields a term p{gh). After a 
substitution h i-^gh we obtain 


|C| 

E Zn.PMG" 




Ck^,qk^ghq^^ 

€2 


g.h^Zcir) fci = l 


As remarked before, every element s € G can be written uniquely in the form 
s = nqi for some i and n € Zc{r). Hence the summation over ki and g can be 
replaced by a summation over s G G. More precisely, we set s = gq/.^. Note 


19 


that Cfcj 
to 


= qk^rqf.^ = srs. With this observation the expression above reduces 


E 


/ > ri232 


{h)F:'3’ 


s^srs,shq.^ ^ ^jj 

42 4 ’ 


h^ZQ^r') s^G 

where the equality follows after a substitution s i—^ and with the help of 

equation (12.51) . This completes the proof. □ 


Analogously to the automorphisms for the abelian models, we now define 
linear maps x/j(A) of 21. Choose again a semi-infinite ribbon ^ and let 
denote the corresponding multiplets, where is the first part of the ribbon, 
consisting of n triangles. Then we set for any local observable A: 


XIM) := lim ^ F/f A (f//)* . (5.3) 

n—¥cci ^ V sTi / 

K 


Note that (assuming for the moment that the expression converges) this defines 
a linear map defined on a dense subset of 21. Since it is bounded, it can be 
extended to 21. This extension will also be denoted by xu- So it remains to be 
shown that the expression indeed converges (in norm, even). The main idea is 
similar as in the abelian case, only here we have to use Lemma l5.II 

Lemma 5.2. Let xu be as above. Then the limit on the right hand side of 
equation (lOD converges. We have the following properties: 

1. for A e %oc, X/j(A) =Y.kHn^ ^ ^^9 enough; 

xij{i) = hji; 

3. X 7 j(A) = 5ijA if A is localised away from the ribbon; 

4- XuiAB) =J2 kXik{A)xkj{B); 

5- Xij{A)* = XJi{A*)- 

Proof. We show the first property. The others then follow straightforwardly 
using orthogonality and completeness for the multiplets, as well as the definition 
of Xu- Consider A G Let N be such that supp(A) D (^„ \ fn) = 0 for all 
n > N. The idea is to decompose the ribbon fn as fn = where ^ = f,n\f.N- 

Let us now write Xiji^) for Y.k (^4^) > = G, 6 = (Cn\?w)- 

By Lemma [5.11 locality, and equations (EH) we get 


x7AA) = E E E F/^^F^^^A(F/^^F^^^)* 

K L M 

= E E E 

K L M 

L M 

= x?AA). 

From this it is clear that the limit in equation (15.31) converges for operators 
A G 2lioc- As mentioned the other properties are easy to verify. □ 
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Note that the properties stated in the Lemma are precisely those that one 
needs to define an amplimorphism. Note that there are n = \C\ dim(p) pairs 
I = We then define a map y : 21 ^ M„(2l) by setting [x{A)]ij = Xij{A). 
It follows that X is an amplimorphism. In addition, it is localised in any cone A 
that contains the ribbon ^ used in the definition of x/j, in the sense that x(^) 
is the matrix with entries A on the diagonal and otherwise zeros, if A G 2l(A'^). 
It remains to show that the amplimorphisms are transportable. Here we show 
that the charges can be transported over a finite region. 

Lemma 5.3. The amplimorphisms x constructed above are transportable over 
finite distances. 

Proof. Fix a semi-infinite ribbon f. We demonstrate how we can move the 
charge at the endpoint of the ribbon around. That is, let ^ be a ribbon, such that 
is again a semi-infinite ribbon. Define V G M„(*8('H)), where n = \C\ dim(p), 
by having entries V/j = Fl'^. From equations (15.11) it follows that V is unitary. 

The claim is that Vx{A)V* = x{A), where the amplimorphism in the right 
hand side is defined with respect to the ribbon f. This can be verified for local 
operators A, by carrying out the matrix multiplication, and using Lemma 15.11 
together with [TJ of Lemma 15.21 □ 

The case of cone transportability is more complicated. Because of the lemma 
it is enough to consider two semi-infinite ribbons and ^2 starting at the same 
site. Let denote the ribbon consisting of the first n triangles. For each n, 
choose a ribbon that goes from the endpoint of ^ 2 ,n to the endpoint of in 
such a way that as n goes to infinity, so does the distance of the ribbon to the 
endpoint of ^ 1 . Now define a unitary operator Vn = n’ where „ is 

the unitary obtained from the multiplet Now if A is local, it follows that 

14iXi(^) = A 2 (A)y„ for all n large enough. This can be seen by the argument 
in the proof of the above lemma. 

This gives a uniformly bounded sequence of operators, since each of them 
is unitary. By the compactness of the unit ball in the weak operator topology, 
there is a subnet that converges to some operator V. Since multiplication on the 
right with a fixed x(^) is weakly continuous, it follows that V intertwines x(^) 
and x(^) • Ths problem remains to show that V is unitary. There are different 
ways that one might achieve this. For example, one could first try to show that x 
and X s-re irreducible, so that Hom(x, x) must be either zero or one-dimensional. 
If V is non-zero, it then follows that one can choose V to be unitary. The other 
option is to realize both representations as the GNS representation of the same 
state. The vector state with an D in the first component (and otherwise zero) 
is a good candidate. By the independence of the ribbon operators on the exact 
choice of ribbon, this leads to the same state in both representations. If one can 
show that this vector is in fact cyclic, the proof is complete. We will leave this 
issue open for now. 

It is natural to look at the amplimorphisms in the ground state represen¬ 
tation, that is, look at (tto ® id) o x- This amounts to applying the ground 
state representation to each matrix element. In fact, that is what we have been 
doing implicitly above. Note that by localization of the amplimorphisms, for ob¬ 
servables outside the localization region, this representation looks like n copies 
of the ground state representation. Hence it is natural to adapt the selection 
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criterion dSD a bit to allow for this case. It turns out that in the end this 
does not really matter, and one can go back from the amplimorphism picture to 
cone-localised endomorphisms (or representations). This is the content of the 
next theorem. 

Theorem 5.4. Suppose that ttq satisfies Haag duality for cones. Let tt be a 
representation of 2t and n be a positive integer such that for some cone A, we 
have 

n • TTo r 2l(A=) ^ TT ( 2l(A‘=), (5.4) 

where n-wg is the direct sum ofn copies of the representation. Then the following 
hold: 

1. There is an amplimorphism x : 21 —>■ M„(2l^“), localised in A, such that 
we have (ttq 0 id) o x = tt. 

2. There is a morphism p : 21 —>■ 2l^“ such that ttq o p = tt when restricted to 

21 . 

Note that the unitary equivalence in the second point can be used to map inter¬ 
twiners between amplimorphisms to intertwiners between morphisms. In par¬ 
ticular, if X is transportable, so is the corresponding morphism, and p can be 
extended to an endomorphism of 2l^“. 

Proof. dH) Let U : Htt 0^=1 be the unitary setting up the equivalence 
and suppose that A € 2t(A). For A s 2t, define x(^) ■= Utt{A)U*. Note 
that x(^) S 1B(0)L;^ ?^o) — Af„(lB(Ho)), so that the matrix elements Xu (^) S 
*8(Ho). Now let B S 2l(A°) and A € 2t(A). Note that from equation (15.411 it 
follows that x(-B) = diag(7ro(i?),... ,7ro(i3)). 

Remark that A and B commute by locality, so that x(A)x(B) = x{^B) = 
x{B)x{A). Writing out the definitions and comparing matrices element wise, 
it follows that ■ng{B)xij{A) = Xij(A)7ro(H). Hence Xiji^) € 7ro(2l(A'=))' = 
ttq (21(A))", by Haag duality. Since the algebra on the right hand side is contained 
in the auxiliary algebra 2l^“, and x acts trivially on 2l(A‘^), it follows that 
X : 21 —>■ Mn (2t^“) is an amplimorphism. 

([2) Since 7ro(2l(A))" is an infinite factor, one can find isometries Vi, i = 
1,... n generating a Cuntz algebra m- That is, they satisfy V*Vj = SijI and 
Write X for the amplimorphism obtained in part (P), and 
Xij{A) for its matrix elements when evaluated in 21. Define a map p : 21 ^ 2l^“ 
by 

n 

p(A) = ^ Vx^J{A)v;. 

Suppose that A, H e 21. A straightforward calculation then shows 
p(A)p(H) = ^ Vx^M)V*VkXkl{B)Vf 

= '£Vx^M)xji{B)vr 

= p{AB). 

Moreover, p(A*) = p(A)* since Xii(A*) = Xji{^)*■ V B € 2l(A'^), then Xij{B) = 
SijB and B G 21(A)', hence p{B) = B and p is localised in A. 
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Next we show that ttq o p is unitarily equivalent to ttq <E> X- To this end, 
identify ® C" with % = 0"^;^ 'Ho- Define a map U ■.'H ^'Ho hy setting 

n 

• © V’n) = y^Vilpi- 

i-l 

Using the properties of the Vi it is easy to check that U preserves the inner 
product (and hence is an isometry and well-defined). Since = I V 

also has dense range, hence U is unitary. □ 

A few remarks are in order at this point. First of all, the condition of Haag 
duality is only used to obtain an amplimorphism from a representation satisfying 
equation (lOl) . Without Haag duality, one can still obtain an morphism p 
such that ttq o p and ttq ® x unitarily equivalent. However, we then have 
little control over the range of the morphism, and we cannot extend it to an 
endomorphism of the auxiliary algebra without any additional information. 

The theorem shows that in principle one can restrict to the study of localised 
and transportable endomorphisms. Nevertheless, it can be very helpful to look 
at amplimorphisms as well. One reason is that it may be easier to construct such 
amplimorphisms explicitly in concrete models, as we have done above. They also 
provide more information on the symmetries of the model. In particular, the 
vector space C" in 'Ho © C” carries a representation of the symmetry algebra, 
through the symmetry transformations of the multiplets . This plays an 
important role in the analysis of the superselection sectors in [52l [58]. We 
expect that the methods used there to study the category of amplimorphisms 
(and hence, by the theorem above, the category of localised endomorphisms). In 
particular, to define fusion and braiding. This leads to the following conjecture: 

Conjecture 5.5. Let G he a finite group. Then the category A of localised and 
transportable endomorphisms of the quantum double model is equivalent to the 
representation category Rep^I?(G) of the quantum double 'T>{G). 

In [46] we obtained through computer algebra the fusion rules of I?(S' 3 ) by 
taking the composition of certain positive maps. These results are consistent 
with the approach here, since the positive maps there are the traces of the 
amplimorphisms we have defined here. The composition of these positive maps 
yields the trace of the amplimorphism xi © X 2 j where the tensor product is 
defined as in [58] . 

To study this conjecture it would be helpful to go from endomorphisms to 
amplimorphisms. An important question is also if there is a “canonical” way to 
obtain an amplimorphism from an endomorphism. Of course, there is always 
the “trivial” way, since and endomorphism p gives rise to an amplimorphism 
p : 21 ^ 21 © C. In general, let Ti,... T„ G 1B(H) be such that T*Tj = diji 
and TfclX = Then it is easy to check that [x(^)]y := T*p{A)Tj 

defines an amplimorphism. Hence to have a meaningful equivalence between 
the category of amplimorphisms and of localised morphisms, we would have to 
impose some additional conditions on the amplimorphisms. One of them is that 
they transform in the right way, as explained above. To find this symmetry one 
could look for an analogue of the Doplicher-Roberts theorem [21] , which gives a 
group symmetry for bosonic/fermion sectors. In general, one could only expect 
a so-called weak Hopf-algebra symmetry [^ . 
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We have outlined here how one could proceed with an analysis of the su¬ 
perselection structure of the non-abelian quantum double model. Although the 
analysis is not complete, we hope that it is a helpful starting point. It should 
be remarked that here we have provided an explicit model whose ground space 
representation should lead to non-abelian charges. This should be contrasted 
with [^, where the existence of such a representation is taken as an as¬ 
sumption. 


6 An application: topological quantum comput¬ 
ing 

There are different reasons why there has been a huge interest in topologically 
ordered systems in recent years. One of the reasons is that they provide examples 
of new phases of matter, that go beyond the Landau paradigm of symmetry 
breaking. This is not just of theoretical interest - these phases really exist in 
nature. There is also experimental evidence for the existence of quasi-particle 
excitations with anyonic excitations, see for example [10]. A good theoretical 
understanding is therefore very welcome. 

Here we focus on another aspect that has sparked the interest of the quantum 
computation community. One of the goals of quantum computation is to use the 
full power of quantum mechanics to solve computationally hard problems, for 
which a computation on a usual, classical computer is infeasible. To illustrate 
this one can think of a simple spin-1/2 quantum system, with Hilbert space 
C^. The dimension of n copies of such a system scales as 2"’, so if one wants to 
simulate the whole Hilbert space of a n-particle system, one quickly runs out of 
memory in a classical computer. The idea behind quantum computing, which 
goes back to Feynman [22] . is to use the laws of quantum mechanics to solve 
computationally complex problems, or even simulate other quantum systems. 

This is not the place for a full-fledged introduction to quantum computing 
(for this Nielsen and Chuang’s [5T] is a good start), but let us summarise the 
main points. The “quantum memory” of a quantum computer is modelled by 
a (generally finite dimensional) Hilbert space, ususally of the form % — (C^)®" 
for some n. A computation then consists of three steps: 

1. Initialise the system in a known state; 

2. Perform a unitary operation on the system to implement the algorithm; 

3. Measure the result (and if necessary, repeat to get statistics). 

This encompasses classical computing. To see this, we can reformulate each 
computational problem in the calculation of a function / : {0,1}" ^ {0,1}”, 
where for simplicity we take / to be injective. Now let |0), |1) be a basis of C^. 
Then we can identify each {xi) G {0,1}” with a basis vector jx) := \xi) ® ® 

\xn)- We can then define a map Uf hy Uf\x) = Uf \ f{x)), which is unitary 
because / is injective. Hence we can initialize the system in a known state |x), 
apply Uf, and measure to learn something about |/(x)). This is where quantum 
mechanics comes in: it allows us to act with Uf on a superposition^ say of the 
form |x), so that we can learn something about all values /(x) in a 

single operation. This is not possible on a classical computer. 
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Although the main idea is quite simple, it turns out to be very difficult 
to implement this in physical systems in a reliable and scalable way. This 
is were topological quantum computing comes in [SIl [SO]. We have already 
briefly discussed the problem of storing quantum information over an extended 
period of time. Here we focus on the computation part, that is, processing this 
quantum information to implement an algorithm. As mentioned above, this 
amounts to acting with a certain unitary operator on the system. In practice 
this could work, for example, by coupling the system to some external magnetic 
field for a period of time, and let it evolve. The problem is that it generally is 
very difficult to exactly perform the unitary that you want, and not something 
slightly different (because the magnetic field is left on too long, for example). 
The idea is therefore to use “topological” operations to implement the necessary 
unitary. A small disturbance should not change the topological property, and 
hence have no effect on the computation. The braiding operation of two anyons 
for example is independent of the path that the anyons take (as long as the 
paths don’t cross). 

To use this idea, the quantum information on which we want to operate has 
to be encoded using anyons. Let us consider a non-abelian anyon p, which for the 
sake of simplicity we assume to be self-dual: p = p. We can then take n copies 
of it, that is, consider p®". The idea is to create these by pulling pairs of them 
from the vacuum, which is possible because p is its own anti-charge. Because 
the anyon is non-abelian, this can actually be done in different ways. More 
precisely, the possible states are described by Hom(/,,p®”), where t is again the 
trivial charge. This vector space can be given the structure of a Hilbert space. 
It is this space that we use to encode the qubits in. It should be noted that in 
general Hom(i,p®) does not have a nice decomposition as the tensor product 
of n copies of some Hilbert space. Nevertheless, one can embed the state space 
of a number of qubits into this space. The dimension of Hom(i,p®") grows 
exponentially in n. 

Step 1 in the quantum computation scheme can now be accomplished by 
pulling charges from the vacuum in a suitable way. To implement unitary oper¬ 
ations on the qubits, we note that there is a natural representation of the braid 
group Bn on Hom(i,p®”): if T G Hom(t, p®”), we can consider 

TT{bi)T := (/ (g) • • • (g) Sp^p (g • • • (g /) o T, 

where bi is the generator of the braid group that swaps the f-th and {i + l)-th 
strand, while on the right hand side, the Ep^p term acts on the f-th and {i + l)-th 
tensor factors. Note that the result is again in Hom(t, p®"), hence this gives a 
unitary rotation on the encoded qubits, and hence allows us to implement (a 
part of) a quantum computation. 

This is not the end of the story, because it is not clear if every unitary on the 
encoded qubits can be obtained in this way. It is enough for 7r(i3„) to generate 
a dense subset of the unitaries. If this is the case for a certain anyon model, it 
is said to be universal, because it means that each quantum algorithm can in 
principle be implemented on it. Kitaev’s quantum double model is universal for 
a wide range of non-abelian groups [131 [H] , but in general this is a rather special 
feature of a model. If a model is not universal, one can supplement the braiding 
operations with other, non-topological operations to be able to implement each 
unitary operation. Even if the model is universal, one still has to find out which 
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combination of braidings one has to do to obtain a certain unitary. Luckily, this 
can be found efficiently using the Solovay-Kitaev theorem m- This is worked 
out explicitly for a simple anyon model, the so called Fibonacci model, in [7]. 

Finally, after the appropriate braidings have performed, it is time to measure 
the outcome. That is, we have to find out which state in Hom(i, p®”) the system 
is in. This is done by fusing some of the anyons again. Recall that the fusion 
rules are of the form p(Sip = -^ppPfe- The integers essentially says in how 
many ways the fusion of two p-anyons can lead to an anyon p^. It are precisely 
these different ways to fuse n anyons to the vacuum that label the states. Hence, 
by fusion some anyons and observing the outcome (which amounts to doing a 
charge measurement), we can obtain statistics on which state the system was 
in. The charge measurement is where it becomes important to have modularity 
of the category: this allows us to distinguish the charges by pulling pairs of 
charges from the vacuum, move one of the anyons around the region for which 
we want to determine the charge, and fuse to the vacuum again and observe if 
there is anything left or not. 

It should be noted that while topological quantum computing has advantages 
with respect to the stability of the operations, there are also drawbacks. For ex¬ 
ample, in practice it may not be so easy to physically move the anyons around, 
especially over larger distances. This may be circumvented by measurement 
based quantum computation. There, the braiding operations are mimicked by 
doing a series of measurements [B] . One could also restrict to more local opera¬ 
tions, although this will mean that the anyon models are no longer universal [3]. 

To conclude, we have seen that there is a rich class of so-called topologically 
ordered states, for which methods of local quantum physics provide useful tools 
to study such systems. Besides the possible applications to quantum computing 
that we have mentioned, there are also very interesting condensed matter and 
mathematical aspect related to such phases. Finally, also on the experimental 
side the field is very active. It would be good to see if the tools of local quantum 
physics can be further employed to advance progress in this multi-disciplinary 
field. 
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